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Abstract 

Based on [T] we introduce and study a mathematical framework for analysis- 
suitable unstructured B-spline spaces. In this setting the parameter domain has 
a manifold structure which allows for the definition of function spaces that have 
a tensor-product structure locally, but not globally. This includes configurations 
such as B-splines over multi-patch domains with extraordinary points, analysis- 
suitable unstructured T-splines, or more general constructions. Within this 
framework, we generalize the concept of dual-compatible B-splines (developed 
for structured T-splines in 0 )- This allows us to prove the key properties that 
are needed for isogeometric analysis, such as linear independence and optimal 
approximation properties for /i-refined meshes. 


1. Introduction 

In isogeometric analysis, as it was introduced in [3], the physical domain is 
usually parametrized smoothly by tensor-product B-splines or NURBS, there¬ 
fore it has to be diffeomorphic to a rectangle. When dealing with geometrically 
complex domains, there is a need for a more flexible geometry representation. 
Several approaches have been proposed so far. We would like to point out sev¬ 
eral common strategies that can be studied within the framework we present 
in this paper. One is to use a multi-patch representation of the domain, i.e. 
to split the domain into rectangular or hexahedral boxes. Starting from the 
framework presented in [1] , there exist several contributions on multi-patch rep¬ 
resentations in isogeometric analysis [3 ig H 0 0 CD]. In most constructions, 
the basis functions are only C° across patch interfaces. We want to point out 
the recent paper m, in which the authors present a construction to enhance 
the smoothness across patch interfaces away from extraordinary features. 

Another way to handle complex geometries is to use a specific set of non¬ 
tensor product function spaces, such as unstructured T-splines (see mm as 
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well as m)- In both cases, the resulting function spaces are very similar to the 
ones proposed in Sections |5.3| and |5.4[ We however follow a slightly different 
approach to represent the parameter domain and define the function spaces over 
it. 

Smooth constructions over general quad meshes can be introduced by means 
of subdivision surfaces. They have been studied in the context of isogeometric 
analysis in [TU [151 IE] ■ We do not go into the details of subdivision based con¬ 
structions, but want to point out the relation to the manifold based framework. 
Function spaces based on subdivision schemes that generalize tensor-product 
B-splines, such as Doo-Sabin (see m) for biquadratics or Catmull-Clark (see 
[18]) for bicubics, can be interpreted as spline manifold spaces. In both cases, 
the function space is spanned by standard tensor-product B-splines as well as 
certain special functions near extraordinary points. The special functions are 
piecewise polynomials over infinitely many rings of quadrilateral elements. 

In this paper we consider spline manifold spaces as a theoretical tool in iso¬ 
geometric analysis. Note that it is not the purpose of this study to introduce a 
new way to define or implement unstructured spline spaces. The aim is rather 
to develop a framework to study the theoretical properties of a wide range of 
spline representations over non-rectangular parameter domains based on mani¬ 
folds. Concerning spline manifolds we refer to [T], where the authors propose an 
abstraction of the concept of geometry parametrizations that is based on a man¬ 
ifold structure. In the proposed framework, the underlying parameter domain 
forms a manifold, thus allowing to represent more general geometric domains. 

In the following we briefly review the concept of a geometric domain which 
is parametrized by spline manifolds, and discuss its main features. We consider 
a domain E C M" which can be interpreted as a d-dimensional manifold with 
d < n. We mostly focus on a two-dimensional planar domain S C K^, a surface 
E C or a three-dimensional volumetric domain E C R^. In an abstract 
setting, a manifold E is defined by an atlas, i.e., a family of charts E^ such that 

N 

i=l 

together with suitable transition maps between intersecting charts E^ and Ej. 
To define E as a spline manifold, we must define a family of open parameter 
subdomains fli, together with a spline space on each subdomain, and each 
chart is the image of a spline parametrization —>■ E^. The union of 

the subdomains Vii forms an unstructured parameter domain 17, and suitable 
transition functions between the subdomains, together with a relation between 
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the corresponding spline spaces, endow 17 with a manifold structure, which 
is inherited by S. Notice that to cover completely the domain S7, the open 
parameter subdomains 17^ overlap. Hence, in order to cover the extraordinary 
points (or edges) it is necessary that some of the subdomains 17^ also have an 
unstructured configuration, and non-tensor product spline functions need to be 
defined. This is in contrast to traditional multi-patch representations, where 
only the patch boundaries intersect, and the (closed) domain is defined as the 
union of the closure of the patches. 

In Section we recall and study the notion of a parameter manifold as 
presented in [T]. We introduce the mesh on the parameter manifold and cor¬ 
responding mesh constraints in Section In Section we introduce spline 
manifold spaces based on tensor product B-splines and relate them to existing 
constructions [101 ns Ha [13]. Finally, we develop the framework of analysis 
suitable spline manifold spaces in Section where we extend the notion of 
dual-compatibility to B-spline manifolds. In the main part of the paper we as¬ 
sume that the manifold has no boundary. For a more detailed study of spline 
manifolds with boundary see Appendix [Appendix A| We conclude the paper 
and present possible extensions in Section 

2. Parameter manifold 

Before we can define unstructured spline spaces on manifolds we need to 
introduce the abstract representation of the parameter domain, the so-called 
parameter manifold. The following definitions are taken from [1] . The definitions 
are valid for arbitrary dimension d, but we will mostly consider d = 1, 2, 3. 

Definition 1 (Proto-manifold). A proto-manifold of dimension d consists of 

• a finite set (named proto-atlas) of charts LOi, that are open 

polytopes oji C line segments for d = 1, polygons for d = 2 or polyhe- 
dra for d = 3; 

• a set of open transition domains such that Wij C uJi and 

■> 

• a set of transition functions that are homeomorphisms 

'ipij : LOij —>■ ujpi fulfilling the cocycle condition o in 

Wij n LOi^k for all i, j,k = 1,..., N; 

• For every i,j, with i ^ j, for every 17 oji and Cj G dwj^i n Wj, 

there are open balls, and Vi^., centered at Ci and Cj, such that no 
point of n uij^i is the image of any point of H ujij by 
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Figure 1: Visualization of the cocycle condition 


In Figure [l] we visualize the cocycle condition, depicting three domains oji, 
LOj and LOk and respective transition functions V’i.ji '4’j,k and '4’i,k- The hatched 
regions represent corresponding transition domains. 

In general, the transition domains coij may be empty. It follows directly from 
the cocycle condition that the transition function is the identity function 
on LOi^i = oJi and that = tpj^i for all i,j. 

Remark 1. The last condition in Definition^is taken from m- It guarantees 
that the proto-manifold actually represents a manifold, i.e. that there are no 
bifurcations of the domain. Note that this condition is not necessary for the 
following definitions. However, if the condition is omitted, the resulting object 
is not a manifold anymore. 

In the next section the charts Wi will be used as (local) parameter domains 
to define splines. By merging and identifying the charts of the proto-manifold 
we obtain a manifold, which will serve as the parameter domain in our setting, 
thus the name parameter manifold. 

Definition 2 (Parameter manifold). Given a proto-manifold, the set 



( 1 ) 


is called a parameter manifold. Here [J denotes the disjoint union, i.e.. 


U = {[C,*], C e = I,...,iV} 


2=1,...,V 


and the equivalence relation ^ is defined for all Ci ^ Cj ^ 


[Ci,*] ~ [Cj, j] = Cj- 
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We denote by T^iiCi) & the equivalence class corresponding to Ci S 

To motivate the definition (and wording) above, we recall that is indeed a 
topological manifold, endowed with the natural structure (see [1] and [H]): 
is a one-to-one correspondence between each uji and Tri(uii) and plays the role 
of a local representation, its inverse 7r“^ is the coordinate chart in the classical 
language of manifolds. In the following we use the notation and 

= Tri(uJij), where fli, C fl. Hence we have H Hj. The 

parameter manifold H can be interpreted as a topological manifold. 

Remark 2. To be •precise, a parameter manifold is a class of piecewise smooth 
manifolds, which is a sub-class of topological manifolds and a super-class of 
smooth manifolds. It is similar to the concept of piecewise linear manifolds (see 
Depending on its local structure, the parameter manifold is either or 
locall'y. In general the parameter manifold is globally and piecewise C°°. 

We assume that there exists a metric on H, that allows us to introduce the 
usual Lebesgue space lf[VL). In Sectionj^we present in more detail the function 
spaces that are necessary for the analysis. 

Having defined a parameter manifold O, we can define a mesh T on H in 
Section and a function space with a piecewise construction over the mesh in 
Section IH 

3. Mesh and mesh constraints on the parameter manifold 

We introduce the general concept of a mesh T on a parameter manifold H in 
Section |3.1| and then develop the specific configurations we consider in Section 
|3.2| Following that, we present some example configurations and study the 
meshing of complex geometries in Section [3.3[ 

3.1. Mesh on a parameter manifold 

First we define a proto-mesh on the charts. 

Definition 3. A proto-mesh 

with Ti = {q C wj (2) 

is a collection of sets, where 

• each set is composed of open polytopes q, called elements', these are 
intervals, quadrilaterals, hexahedra, etc., depending on the dimension d = 


5 


1, 2,3,..respectively, and each set is a mesh on Wi, i.e. the elements 
are disjoint and the union of the closures of the elements is the closure of 

• for every i,j, ujij is the interior of the union of the closure of elements of 
Tf, and 

• the transition functions map elements onto elements, i.e. 

yq & Ti, G Tj. (3) 

Furthermore, we have the following. 

Assumption 4. The transition functions are piecewise d-linear mappings with 
respect to the mesh. 

The proto-mesh naturally defines a mesh on fl. 

Definition 5 (Mesh on fl). We define the mesh on the parameter manifold fl 
as 

T={Qcn: Q = Tri{q),qGTi,i=l,...,N}. (4) 

Remark 3. Thanks to ([^, the set 7' in ^ is indeed a well defined mesh on 
n. The elements of T are subsets of Tl that fulfill the standard properties of a 
mesh, i. e. the elements are disjoint, the union of the closures of the elements is 
the closure offl. Since Q is a topological manifold the notion of the closure of 
elements, boundary edges, faces, etc. is well defined and derives from the local 
definition on each chart and the equivalence relation given by the transition 
functions. 

3.2. Structured and unstructured meshes on the parameter manifold 

In the following we classify and restrict to specific but relevant charts, de¬ 
pending on their local mesh topology. 

We will say that several objects share a common object A, if A is contained 
in the closure of each of the objects. Mesh objects and properties defined on 
the charts can be carried over to the mesh T on fl. We will not distinguish 
between the mesh T on and the proto-mesh defined on the charts, unless it is 
necessary. The relations between vertices, edges, faces and elements are stated 
within the global mesh. Moreover, these geometric objects are always assumed 
to be open. 

Definition 6 (Structured chart). A d-dimensional chart uji is called a structured 
chart if 


6 


(A) uji is a d-box, that is 


d 

uJi = 

i=i 

and Ti is a box mesh (see, e.g., EH); 

(B) for every j, ujij is a d-box or the union of d-boxes with disjoint closures; 
and 

(C) for every j, if ujj is a structured chart then the transition function tpij : 
LOij —>■ ujj^i is an affine mapping (a linear polynomial) in each connected 
component of ujij. 

Considering unstructured charts we need to distinguish three different types, 
depending on the dimension and on the topological structure. 


Definition 7 (Unstructured vertex chart, two-dimensional). A two-dimensional 
chart LOi is called an unstructured vertex chart corresponding to an extraordinary 
vertex of valence ^ 4 if 


(A) oji is formed by a ring of ki conforming quadrangular segments Si^£, with 

i = around where each segment has a mesh ai^i which is 

topologically equivalent to a box mesh, and the mesh ti is given as the 
union of the meshes ai^i on the segments Si_^; 

(B) LOi, except for the extraordinary vertex is covered by the transition 
domains LOij with structured charts, i.e. 

U 

j = l,...,N 

ujj structured 

and the transition domains tOij are given by the union of segments; and 

(C) if LOj is an unstructured chart with i ^ j, then tOi j = 0. 


To give more insight into the definition of an unstructured vertex chart via 
its segments, we present some example configurations, which are all constructed 
from the same three conforming segments. Figure [2(a)| depicts a mesh where 
each segment is covered by a single element. Figure [2(b) | depicts a conforming 
mesh and Figure 2(c) contains two hanging vertices, hence it is a non-conforming 
mesh. This figure also tells that the non-conformity of the mesh over the un¬ 
structured vertex chart may have two reasons, either the meshes on two neigh¬ 
bouring segments do not match (upper right hanging vertex) or the T-node is 
already present within the mesh ai^£ on the segment (lower left hanging vertex). 
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(a) 


(b) 


(c) 


Figure 2: Examples of meshes on a ring of 3 conforming segments around a two-dimensional 
extraordinary vertex of valence 3 

Definition 8 (Unstructured edge chart). A three-dimensional chart oJi is called 
an unstructured edge chart of valence ki ^ A li 

(A) uii = (hi X u)i, where the two-dimensional chart hi is a two-dimensional 
unstructured vertex chart, as in Definition!^ with an extraordinary point 
^i of valence ki, partitioned into two-dimensional segments Si^( for i = 
1,..., /ci; uii is an interval ioi = ]ai_ 3 , bi^sl; each three-dimensional segment 
Si,e = Si^i X uji has a mesh which is equivalent to a three-dimensional 
box mesh, and the mesh is again given as the union of all meshes ai^t 
on the segments; 

O 

(B) Wi, except for the extraordinary line x Qi, is covered by transition do¬ 
mains Wij with structured charts, i.e. 


O 




LJj structured 


and the transition domains (jOij are given as the tensor-product of the 
union of two-dimensional segments with an interval in the third direction; 

(C) for structured charts ijjj the transition function ipi,j is linear in the third 
coordinate Cs; and 

(D) if LOj is an unstructured edge chart with i ^ j and Uij ^ 0, then there 
exists a structured chart Wfc such that iVij C 


Definition 9 (Unstructured vertex chart, three-dimensional). A three-dim¬ 
ensional chart OJi is called an unstructured vertex chart corresponding to an 
extraordinary vertex of valence ki if 












(A) LOi is formed by a ring of ki conforming hexahedral segments with 

£ = around where each segment has a mesh Ui^i which is 

topologically equivalent to a box mesh, and the mesh is given as the 
union of the meshes on the segments 

(B) LOi, except for the extraordinary vertex is covered by the transition 
domains Wij with unstrnctured edge charts, i.e., 

U 

3 = 1,...,N 

cjj unstructured edge chart 

and the transition domains are given by the union of segments; and 

(C) if LOj is an unstructured vertex chart with i ^ j, then lotj = 0. 

Note that from Definition |^(B) and Definition]^ (B) it follows that for every 
unstructured vertex chart in 3D the transition domains with structured charts 
cover everything, except for the extraordinary features (union of extraordinary 
vertex and extraordinary edges). In particular, the interior ef every element is 
covered. 

We would like to point out that Definitions [!l(c).§ (D) and (C) are not 
strictly necessary, but technicalities to simplify the definition of spline manifold 
spaces in Section They guarantee that unstructured vertex charts do not 
overlap and that unstructured edge charts only overlap in the vicinity of an 
unstructured vertex. Notice that it is not possible that one element of the mesh 
T is adjacent to two unstructured vertices. In the context of meshing, this 
is not a severe restriction as we point out in the following section. Moreover, 
each unstructured vertex and each unstructured edge corresponds to exactly 
one unstructured vertex chart or unstructured edge chart, respectively. 

Assumption 10. We assume that each chart tOi can be either structured as in 
Definition or unstructured as in Definitions [7)[8]and[^ 

In the following section we give some example configurations and study in 
more detail the meshing of complex geometries. 

3.3. Example configurations and meshing of complex geometries 

For simplicity, in the classification above we have restricted ourselves to a 
limited number of types of charts, which nevertheless contain most mesh con¬ 
figurations of practical interest. We first give a summary of the different types 
of vertices that can occur (depending on the dimension), which motivates our 
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chart classification. Each unstructured chart is formally defined in such a way, 
that it can be used to cover a certain type of unstructured vertex. Following 
the discussion of the types of vertices, we present some simple example con¬ 
figurations and discuss the issues of meshing with manifolds. Without beeing 
thorough, we present configurations that can be covered and such that cannot. 
Since we do not want to go into the details of meshing we refer to the following 
literature for quad-meshing [22] and hex-meshing |22l [9] [23] . 

For ID domains (intervals, planar or spatial curves) there are only regular 
vertices, which are equivalent to the knots in the classical B-spline language. 
Hence no unstructured charts are needed. 

For 2D domains we consider three types of vertices. A 2D vertex can be a 
regular vertex, a hanging vertex (a T-node), or an extraordinary vertex (see Ta- 
blej^and Figure]^. Both regular and hanging vertices are covered by structured 
charts, although they may also belong to an unstructured chart. Unstructured 
vertex charts cover a neighborhood of an extraordinary vertex. 


structured vertex regular vertex 

(Figure 

3(a) 

hanging vertex 

(Figure 

3(b) 

unstructured vertex extraordinary vertex (Figures 

3(c) 

and 

3(d) 


Table 1: Classification of vertices in 2D 



(a) (b) (c) (d) 


Figure 3: Different types of vertices in 2D ((a) regular, (b) hanging, (c) extraordinary valence 
3, (d) extraordinary valence 5) 

For 3D domains we have three types of edges and four types of vertices. 
A 3D edge can be a regular edge, a hanging edge, or an extraordinary edge of 
valence fc 7 ^ 4. A 3D vertex is called a regular vertex, if it is shared by regular 
edges only; or a hanging vertex, if it is shared by regular and hanging edges. 
Otherwise, the vertex is called an unstructured vertex. The notion of unstruc¬ 
tured vertices contains both vertices of extraordinary edges, so called partially 
unstructured vertices, as well as fully unstructured vertices. To be precise, a 
partially unstructured vertex is a vertex that is shared by exactly two extraor- 
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dinary edges that can be covered by a single unstructured edge chart. If this 
is not possible, the vertex is a fully unstructured vertex. For a complete list of 
types of vertices and edges in 3D see Table Regular and hanging edges as 
well as regular and hanging vertices are covered by structured charts. Hence, 
we have considered two classes of three-dimensional unstructured charts: un¬ 
structured vertex charts (covering fully unstructured vertices) and unstructured 
edge charts (covering partially unstructured vertices and extraordinary edges). 
Unstructured edge charts are tensor products of an unstructured vertex chart 
in two dimensions and an interval in the third dimension: there is an inner 
sequence of extraordinary edges and all other interior edges are regular or hang¬ 
ing. Similar to the two-dimensional case, an unstructured vertex chart in 3D 
covers a neighborhood of the fully unstructured vertex. 


structured edge 

unstructured edge 
structured vertex 

unstructured vertex 


regular edge 

hanging edge 

extraordinary edge 

regular vertex 

hanging vertex 

partially unstructured vertex 

fully unstructured vertex 


Table 2: Classification of edges and vertices in 3D 


Let us consider an unstructured vertex chart uii composed of ki segments, 
where each segment is meshed with only one element. From Definition it 
follows that for any structured chart ujj with uJij ^ 0 the transition domain 
ipijiujij) = ujj^i is a box-mesh. The same holds for tOij, which is then formed 
by one or at most two quadrilateral elements. There must be (at least) ki subsets 
ujij, each one formed by two adjacent elements, in order to cover uji \ 

Figure [4(a) I depicts a two-dimensional unstrctured vertex chart uJi of valence 

3 and a structured chart, as well as the corresponding transition domains. Since 
the transition domains are open, it can be observed easily, that in this case 3 
structured charts are needed to cover uJi. Figure [4(b)] depicts a three-dimensional 
unstructured vertex chart uJi of valence 4 as well as an unstructured edge chart of 
valence 3 and corresponding transition domains. Again, we may observe that uJi 
(except for the extraordinary vertex) can be covered by the transition domains of 

4 unstructured edge charts. Moreover, (except for the extraordinary vertex 
and edges) can be covered with 6 structured charts. Figure 4(c) depicts an 
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Figure 4: Different types of unstructured charts 


unstructured edge chart of valence 5 and two possible transition domains. Note 
that all the transition domains are mapped box-meshes, due to the assumptions 
above. In all the examples presented here we assume that for each unstructured 
vertex chart each segment is meshed with exactly one element. This is not 
necessarily the case, as we presented in Figure Note that, by definition, the 
transition domains with structured charts can cover no more than 2 segments 
of an unstructured chart. 

Note that the types of charts we consider are sufficient to represent most 
meshes of practical interest. Given an arbitrary quad- or hex-mesh without 
hanging vertices or edges, a global bisection of the mesh can be covered by 
structured charts, as in Definition and unstructured charts, as in Definitions 
and 1^ 

This statement becomes clear, when looking at the types of vertices and 
edges that can occur in the bisected mesh. We consider only 3D meshes in the 
following. We show that every vertex of the bisected mesh can be covered by a 
valid chart of one of the three categories. Every vertex of the initial hex-mesh is 
one of the three: a structured vertex, a partially unstructured vertex or a fully 
unstructured vertex. The bisection of the initial mesh introduces new vertices 
from midpoints of edges, faces and hex-elements (one point for each edge, face 
and element). The midpoints of structured edges, faces and elements become 
structured vertices. The midpoints of unstructured edges become partially un¬ 
structured vertices. Hence no new fully unstructured vertices are introduced. 
All structured vertices can be covered by structured charts, all partially un¬ 
structured vertices can be covered by unstructured edge charts and all fully 
unstructured vertices can be covered by unstructured vertex charts. It is easy 
to see that, due to the bisection, the closure of any element can contain at most 
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one extraordinary vertex. All other assumptions are trivially fulfilled. 


4. Spline manifold space on the parameter manifold 

In this section we introduce spline spaces over the mesh T on fl, in short 
spline manifold spaces. 


4-.1. General spline manifold spaces 

Again based on [T], we define a spline space on a parameter manifold using 
the charts. This is achieved by defining proto-basis functions on the proto-mesh 
(Definition 11) and transfering them onto the parameter manifold using the 
equivalence relation induced by the transition functions (Definition |13[ ) . Each 
chart (jJi plays the role of a local parameter domain. 

Definition 11 (Proto-basis functions). We define a proto-basis as a set 


{{bAi ■ COi -A R}AiGA}i=l,...,Ar ) (5) 

where all proto-basis functions bAi, with S At, are linearly independent 
functions dehned on We further assume that for each i = 1,...,A and 
Ai G Ai the function bAi fulfills 


lim 6 a, (0 = 0. (6) 

Assumption 12. The proto-basis functions bAi with Ai G Ai are piecewise 
polynomials with respect to Tj, i.e. ^A^lg is polynomial for all q G Ti- 

Definition 13 (Spline manifold space). For each i = 1,..., A and Ai G Ai we 
define i?Ai : D —>■ K such that 


^AilOi = bAiOTrA 

J^Aila\Q, = 0, 

and set 

Bi = {Ba- : fl —>■ R, Ai G Ai}. 
Furthermore, we introduce the global index set 


(7) 

( 8 ) 



where the equivalence relation « is defined as follows: given [A.,*] and [AjAI 
in lJi=i N-^i, then [Ai,i] « [Aj,j] if and only if the two functions i?Ai S Bi 


13 



and Ba- G Bj coincide. Therefore, Ba : —?> M is well defined for A G and 

we set 

B = {Ba,A€A}= IJ B,. (10) 




Moreover, we set 


Bi — {^Aloi) A G A-i}, 


where 


Ai = {A G A : supp(i?A) n supp(i?Aj ^ 0 for some A^ G Ai}. 


( 11 ) 


( 12 ) 


The set Bi is the restriction of B onto containing all the functions in Bi and 
also the restriction of any function whose support intersects supp(yBi) but is not 
included in Bi. The functions in Bi can be pulled back to the chart 

Finally, the span of functions in (10) is the spline manifold space, the spline 
space on the parameter manifold fl, denoted by 

5 = spanjSA, A G A} . (13) 

With some abuse of notation, for the global index A G A we will say that 
A G Ai if there exists an index Ai G Ai such that its equivalence class through 
« is equal to A. 

To distinguish different types of functions, we introduce the notation 


As — Ai, 

ijJi is structured 


— 


A.q, - 


U A, I \ As 

yuji is unstr. edge 


(14) 


u 


A, \ (As U Ae) , 


\<jJi is unstr. vertex 


where As is the index set of structured functions, Ae is the index set of edge 
(or partially) unstructured functions and Av is the index set of vertex (or fully) 
unstructured functions. It is clear that 


As n Ae = As n A« = Ae n Ai, = 


(15) 


Remark 4. Equation together with ([^ guarantees that the functions are 
globally continuous. If discontinuous functions are allowed, the condition (§ 
can be omitted. 
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4-2. B-Spline manifold spaces 

We assume that each function that is completely supported in one chart is 
a function of that chart. Moreover, we assume that each structured chart is 
completely covered by the supports of its functions. 

Assumption 14. For each i = 1,..., N and each A G A we have that if 

supp(i3A) ^ 

then A G Ai. If Wi is a structured chart as in Definition then 

y supp(Ba,) = Di. 

AiGAi 

Moreover, if the support of a function is structured, then there exists a 
structured chart that covers the support. 

Hence we conclude that the set Bi contains all functions that are completely 
supported in and the set Bi contains all functions that have a support inter¬ 
secting with rii. 

For simplicity we assume that the functions have the same degree p in each 
direction on each structured chart. 

Assumption 15. If oji is a structured chart as in Definition|^ then each func¬ 
tion Ha in Bi is a tensor-product B-spline of degree p when restricted to uii, 
i.e., there exist (local) knot vectors Sa.i.i, ..., SA,i,d such that 

HA0 7i’i(C) = &[SA,j,i](Ci)--A[SA,*,d](Cd), VAGAi,VCGWi, (16) 

where 6 [S] (C) is the univariate B-spline with local knot vector S. In this notation 
the degree p of the B-spline is given implicitly by the length of the local knot 
vector. 

Finally, the previous assumption is extended to unstructured edge charts, 
since they behave like structured charts along the third parametric direction, 
and like unstructured vertex charts along the first two parametric directions. 

Assumption 16. If is an unstructured edge chart as in Definition then 
each function Ha in Bi is a product of a bivariate function and a B-spline of 
degree p in (^3 when restricted to uii, i.e. there exist a function /3 a, i and a local 
knot vector such that 

-Sa 0 7ri(C) =/3A,*(Ci.C2)&[5A,i,3](C3), VA G .4*,VC G Wj. (17) 
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Remark 5. From the definition of the unstructured charts ujj, and under As¬ 
sumption^^ we have that the functions Ba, with A G Aj, are fully defined by 
their restrictions to the transition domains ujj^i with structured charts. On each 
ujj^i the functions inBjCiBi are equivalent to (mapped) tensor-product B-splines 
according to (161. 


4-3. Isogeometric function spaces using spline manifolds 

We consider a domain E C M" which can be interpreted as a d-dimensional 
manifold with d < n. The most interesting cases, and the ones we focus on, 
are a two-dimensional planar domain E C a surface E C or a three- 
dimensional volumetric domain E C 

Definition 17. The physical domain E is given by a spline manifold parametriza- 
tion G G (5)", with 

G : n ^ E. (18) 

In practical situations, this parametrization is defined by associating a con¬ 
trol point to each function in B. Notice that the geometry E inherits the man¬ 
ifold structure of the parameter manifold fl. Indeed, the parametrization G 
can be considered as a piecewise defined function, where G^ : 17^ —>■ E^ is a 
tensor-product B-spline parametrization for every structured chart 17^. Then 
we have 

N 

S = U^- 

i=l 

where the charts E^ form an atlas of E. 

Then, the isogeometric function space over the manifold E is given as follows. 

Definition 18. Given a spline manifold parametrization G G (5)" as in Defi¬ 
nition [m we define on the manifold E the isogeometric function space as 

V = {/ : E —>• K, such that / = / o G“^ for / G 5}. (19) 

The isogeometric space is well-defined, if the geometry parametrization is in¬ 
vertible. We could move here the assumption on detVG^VG. 

Similar to the spline manifold space S itself, the isogeometric space V can 
be interpreted as a piecewise defined function space. Indeed, each function in 
Bi can be composed with G“^ to define the corresponding function in E, with 
its support contained in the chart E^. 
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4-4- Relation to existing constructions 

We note that several constructions of unstructured spline spaces existing in 
the literature fit in the framework of B-spline manifolds, in some cases with 
minor modifications. For instance, the definition of unstructured T-splines as 
presented in mm for quad-meshes and hex-meshes, the G^-continuous un¬ 
structured T-splines as presented in m. as well as multi-patch B-splines with 
enhanced smoothness as developed in m- The development of an abstract 
framework for the construction of these spaces can serve as the starting point 
for a deeper mathematical analysis of their properties. In the following we detail 
how these three particular examples fit into the framework of B-spline manifolds. 
In all three cases, the idea is to split the set of basis functions into structured 
and unstructured basis functions, and introduce a set of charts covering the 
whole mesh. In this context, a function is called structured, if its support is 
covered by a structured mesh (possibly with hanging nodes). Otherwise it is 
called unstructured. 


T-splines over unstructured meshes m- In this configuration the mesh T we 
consider on is the Bezier mesh and not the T-mesh. Then, we can define 
the set of charts by taking the support of each structured function as a single 
structured chart, while the support of each unstructured function is taken as 
one unstructured vertex chart. In the construction of m it may happen that 
one element contains two extraordinary vertices. In this case two unstructured 
vertex charts overlap, violating Definition]^ (C). This condition is a technicality 
which simplifies the mathematical framework, and could be removed. It can 
also be fulfilled with one level of refinement of the T-mesh, as we explained 
in Section |3.3| In a similar way, the constructions for trivariate functions as 
presented in m fit, with some technical restrictions, into the present framework. 


-smooth T-splines over unstructured meshes m- This construction is similar 
to the construction in m- The support of each structured function can again 
be interpreted as a single structured chart. For each extraordinary vertex we 
can define one sufficiently large unstructured chart u>j, such that it covers the 
support of all functions that are non-zero at the extraordinary vertex. In this 
case, the unstructured functions do not fulfill Assumption |15[ since their degree 
is increased in the vicinity of the extraordinary vertex, and their restriction to a 
structured chart is not a B-spline basis function as in (16), but a suitable linear 
combination of B-splines of higher degree. 
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Multi-patch B-splines with enhanced smoothness m- The authors propose a 
construction based on a multi-patch representation of the domain. This can be 
interpreted as a manifold structure with closed charts (patches), that intersect 
only at the boundary. However, since there is a one-to-one correspondence of 
parameter directions along any shared boundary between two patches, the multi¬ 
patch representation can be transformed naturally to a parameter manifold 
representation by simply enlarging the patches along the parameter direction 
crossing the boundary to obtain open charts. Note that in general it is necessary 
to define more than one chart to cover one patch. Again, there are unstructured 
functions at the extraordinary vertices, and unstructured vertex charts have to 
be introduced to cover their support. 


5. Analysis-suitable spline manifold spaces 


We introduce in this section the conditions for the construction of analysis- 
suitable B-spline spaces on manifolds, that is, spaces that have good properties 
for the solution of differential problems. The key tool for this construction is 
the definition of a (stable) dual basis, which is a set of functionals {Aa, A € A} 
such that 

Aa{Bjs_i) = Saa', VA, A'gA, 

where ^aa' represents the Kronecker delta. A condition for the construction 
of a dual basis for structured T-splines was given in [24l [2] , under the name of 
dual-compatibility. 

We present the construction of dual functionals on the parameter manifold Q 
in Section[5.1[ starting from a proto-dual basis on each chart, and then following 


the scheme of spline manifold spaces introduced in Section [4~T| To guarantee 
that these dual functionals form a dual basis we need to add some conditions 
to the spline manifold. In Section |5.2| we give a dual-compatibility condition 
for spline manifolds, which generalizes the condition in mis] to the unstruc¬ 
tured setting. In this configuration, a global dual basis can be derived from 


the proto-dual bases defined on the charts. Then, in Section 5.3 we present 
an explicit configuration, with a specific construction of the proto-basis func¬ 
tions and the proto-dual basis on unstructured vertex and edge charts. In this 
configuration, which is only C^-continuous at extraordinary features, the dual 
functionals in H can be derived from the proto-dual basis without any further 
modification, which also guarantees the stability of the dual functionals. Finally, 
in Section[5.4|we show the typical application of a dual basis: we prove otpimal 
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approximation properties of the isogeometric space on a simple but interesting 
example configuration. 


5.1. The dual basis for spline manifold spaces 
In the following we use this well-known result. 

Proposition 1. Given a set of finitely many functions {6q,} that are linearly 
independent, there exist functionals Aq. : t K that are dual to the functions, 

i.e., Aq.(6q,') ■ 


As for the definition of the basis functions for spline manifold spaces intro¬ 
duced in Section 4A the starting point is a set of dual-functionals on each chart, 
that exist thanks to the previous proposition and Definition [g 

Definition 19 (Proto-dual basis). We define a proto-dual basis as a set 
{{Aa; : L^{oJi) —t K}AiG.4i}i=l,...,Afi 


where the functionals Aa^ form a dual basis for the proto-basis functions in the 
chart, that is 

AAi(&A') = (^AiA') VAi,A(€Ali. 

Note that the existence of a proto-dual basis is guaranteed by Proposition 
since the proto-basis functions bj^. are linearly independent. We will assume 
that the proto-dual functionals satisfy the following. 


Assumption 20. For any indices [Ai,i] « [Aj,j] that belong to the same 
equivalence class A S A, defined as in §, it holds that Aa; {4> o T^i) = Aaj {4> o 
TTj), V((>eL2(fI). 

The assumption that two equivalent indices are associated to the same dual 
functional is natural, since by ^ they are also associated to the same basis 
function, and it allows the following global definition. 


Definition 21 (Manifold functionals). For any A G A we define the functional 
Aa(((') = AAi((('0 7rj), V(^GL^(D), (20) 


where [Ai,i] is one instance of the equivalence class A. 


From (20), the support of the dual functional A a is contained in 12^. The 
most used dual bases for splines, such as the one by Schumaker [5S] and the ones 
by Lee, Lyche and Mprken [5^ satisfy the stronger condition that the support 
of A A is the same of the corresponding function Ba- 
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Note that in general the dual functionals {Aa, A S A} do not form a dual 
basis for B = {Ba,A G A}. In the following section we introduce a new 
condition that ensures that a dual basis can be derived by a modification of the 
dual functionals Aa- 


5.2. Dual-compatible spline manifold spaces 

For the definition of the dual-compatibility condition we follow and extend 
the recent review paper [27j. Two knot vectors S' = ... ,^p_|_ 2 } and S" = 

... ,^^'+ 2 } overlap if there exists a knot vector S = {^ 1 ,... and two 
integers k' and k" such that and for i = 1 ,... ,p -f 2. 

We generalize the dual-compatibility condition to spline manifold spaces in the 
following definition. 


Definition 22 (Dual-compatible spline manifold spaces). Under the previously 
stated assumptions, the set B defined in (101 is dual-compatible if the following 
conditions hold 

1. for all A G Ai C As and VA' G Ai there exists an index I G {1,... ,c?}, 
such that the knot vectors Sa,*,/ and 'B.A'.i.i as in (16) are different and 
overlap; 

2. for all A G AiDAe, with tOi being an unstructured edge chart, and for all 
A' G Ai, either the knot vectors SA,i ,3 and SA',i ,3 as in 0 overlap and 
are different, or /3a, i ^ /3A',i; 

3. for all A G Ai H Ae, with Wi being an unstructured edge chart, and for all 
A' G (Ai \ Ai) n (Ae U Av) the knot vectors SA,i ,3 and SA',i ,3 as in dli 
overlap and are different. 


Here As, Ae and Ay are defined as in (14). 

m ( as well as the previous papers mw the authors deal only with 
the structured case, and actually the dual-compatibility condition in m cor¬ 
responds to point 1 in Definition 22 The new definition extends the dual¬ 


compatibility condition to an unstructured configuration. As in m, Definition 
[22| gives a sufficient condition for the existence of a dual basis. For that, we 
use two technical ingredients: the first is a univariate L^-stable dual functional 
A[S], such that, if S and S' are overlapping, 


mm']) 


lifS^S', 

OifS^S'; 


( 21 ) 


see, for example, |211 [2^]; the second is the existence of a proto-dual basis on 
each chart, guaranteed by PropositionOur main result follows. 
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Theorem 1. If the set B is dual-compatible, then there exists a set of functionals 
B* = {Aa : —>■ K such that A S A] , 

that is dual to B, i.e., 


VA, A' G A, = Saa'- 


( 22 ) 


Proof. First, we construct dual functionals for the different types of charts. Let 
uJi be a structured chart. Given A £ Ai C As we set 


Aa = A[SA,i,i] G ... G A[sA,i.d], 


(23) 


where the knot vectors Sa.i,* are given as in (16). Let uJi be an unstructured 
edge chart. Given A G AiCi Ae, let A^.a C Ai n Ae be the set of indices A' 
such that SA,i ,3 = SA',i ,3 as in Since the set Bi is linearly independent, 

the same holds for the set 


{/3a'.* : A' G Ai,A} (24) 

and by Propositionthere exist functionals {Aa'.i : A' G A^^a} that are dual 
to (24). We define the proto-dual functional Aa : Wi —>■ K 

Aa = Aa.i G A[SA,i,3]) 


Finally, let uJi be an unstructured vertex chart. Gonsider Ai H Av, and again by 
Proposition [T] the set Bi, which is assumed to be linearly independent, admits 
a proto-dual basis {Aa : A G Ai n A.^}. Having defined A a for all A G A, we 
construct the corresponding Aa as in ( p0| ). 

It is easy to check that, by the construction above and Definitionj^ it holds 


VAg A.,VA'g (A.UAeUA.), Aa(Ha') = ^AA', (25) 

VAg Ae,VA'G (AeUA,), Aa{Ba') = 6aa', (26) 

VAg A„,VA'g A., AAiBA') = Saa'■ (27) 

As already noted, from Definition 21 we also infer 

VA G A,VA'G A\A„ Aa(Sa') = 0. (28) 


In general, the set of functionals {Aa}a6 A is not dual to the set of functions 
B = {Ba}agAj but we can easily fix it by defining, for all if G L^(H), the 
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functionals: 


VAsA, Aa((/)) =Aa (</>), 

VA s A, Aa(</')= Aa(^!>) — ^ Aa(-Ba")^a"(</>), 

A"6(Ai\-4i)nA<, (29) 

VA s Aa(</')= Aa(^!>) — ^ Aa(-Ba")^A"(</')• 

A"GAi\Ai 


By going through all combinations for A, A' S As,Ae,Av and using (25)-(27l 
as well as the condition (28), it follows that the set {Aa}agA as defined in (29) 
fulfills (12^. □ 


5.3. Definition of a basis and corresponding dual basis for a simplified configu¬ 
ration 

The dual compatibility condition of the previous section ensures the exis¬ 
tence of a dual basis, which implies the linear independence of the basis functions 
in B. We now focus on a simplified configuration of B-spline manifold spaces, 
which allows for exactly one extraordinary function in every unstructured ver¬ 
tex chart. We explicitly construct a dual basis within the abstract framework 
of Section |5.2| However, taking advantage of the specific case considered here, 
we obtain dual functionals whose support is the same as the support of the 
corresponding functions. 

In Assumption |23| we restrict our studies to more simple unstructured vertex 
charts, where each segment corresponds to an element of the mesh. Moreover, 
we give an explicit representation for the extraordinary vertex functions. Note 
that one may also define a collection of extraordinary vertex functions that are 
discontinuous at all or some of the element boundaries within the unstructured 
vertex chart. For simplicity, we consider only one continuous extraordinary 
vertex function for each unstructured vertex chart. 

Assumption 23. For each unstructured vertex chart uji, we assume that the 
mesh Ti is the union of the meshes at/ on the segments Si^i, where each ai^i 
contains a single element qi i which covers the whole segment. Moreover, we 
assume that there exists exactly one index A S Mi and the function Ha € Bi 
is named unstructured vertex function. The extraordinary vertex function is 
continuous, has the value 1 at the extraordinary vertex, and for each element 
Qi,e G D there exists a d-linear mapping g ]0, such that 

Ba o irfiC) = bdO ipi^eiC), for all C G (30) 
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where bd is defined as in 


bd{Ci,...Xd) = {i-cir...{i-Cdr, 

= ^[^o](Cl) • ■ ■ ^[“o](C<i)j V(Cl, ■ ■ ■ , Cd) G ]0i j (31) 

with “0 = [ 0 ,... , 0 , 1 ]. 

For each unstructured edge chart uii = uji x Cui, the mesh of the underlying 
bivariate chart Wi is again given such that each segment of the bivariate 
chart is covered by a single element qi^i, the one-dimensional structured chart 
uji is partitioned into a mesh fi = {qi ,..., gmi}) and the mesh on uji is given 
via 

n = {qi,t X qk : (ji^e € h, qu G n). 

Moreover, we assume that for all A', A" G Ai the unstructured edge functions 
are given as in ( [T7| , where /Sa'.i = /3A",i are equal to the same 
unstructured vertex function of dimension d = 2. Hence, for each segment 
qi,i X ]ai^ 3 , bi^Ai there exists a bilinear mapping ipi^i : qi^g —>■ ] 0 , if, such that 

Ba o 7 ri(C) = (&2 o^^i,^(Ci,C 2 )) ^[SA.j, 3 ](C 3 ), for all C G q^^e x ]a*, 3 , 6 ,, 3 [, 

corresponding to the representation in 0 - 

Given an unstructured vertex chart uJi this means that, denoting by j{q) a 
chart index such that is a structured chart and q C there exist knot 

vectors ■ ■ ■, ^A,q,j{q),d of the kind 

[£^_^,nor[r,|V^] (32) 

P+1 times p+1 times 

such that 


Ba °'^j(C) ^[^A,(3,j(g),l](Cl) ■ . • . ■ ^[“A,gj(g),d](Cd)i 

for all C such that '4’j{,q),i{C) G q. We can now define a dual basis explicitly, using 
this representation. A similar representation can be derived for unstructured 
edge charts. 

In the following we define proto-dual functionals Aa and corresponding man¬ 
ifold functionals A a as in Definition which fulfill 

^Ai{BA'f) — ^AiA'-j VAi G Ai, and VA( G Ai- (33) 

Hence they form a dual basis for Ha, with A. G A. For each structured function 
(a function completely supported by a structured chart), the dual functional is 
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given as the tensor-product of univariate functionals (211. For an unstructured 
vertex function we construct a dual functional by adding the contributions of 
each segment of the unstructured vertex chart, each being a structured subdo¬ 
main. For unstructured edge functions, we construct a dual functional anal¬ 
ogously to the basis function in that is, the tensor-product of the dual 

functional of a bivariate unstructured vertex function with the dual functional 
of a univariate B-spline in the third direction. 

Theorem 2. Given a dual-compatible spline manifold spaee let 

B* = {Aa,A G A} 

be a set of functionals Aa : —>■ K such that: 

• if uJi is a structured chart, then for all A G Ai, we have 

Aa((/>) = (A[SA,i.i] 0 ... 0 A[SA.i,d])(</>|ni o TTi), V(/)GL^(n), (34) 


where the knot vectors o,re given as in (16); 

• if UJi is an unstructured vertex chart and Bj^, for A G Ai, is the corre¬ 
sponding unstructured vertex funetion, then let Qi^i = and 

ki 


= V(/. G L^n), 


(35) 


t=i 


where 


^a,iA4>) = (A[5o] ® ...® A[So])((/> o TTj o AiA 
for all (j) G Lf{QiA, with Sq = [0, •. •, 0,1]; and 

• if UJi is an unstructured edge chart and B_x, for A G Ai, is a corresponding 
unstructured edge function, then let Qi^e = T:i{qi^i x ]ai^ 3 , foi^a)), and 

ki 


'^a('^) = E Aa.z/(<(>|q„j, V(^ G LA^), 


(36) 


‘ e=i 

where 

^a,iAA = (A[So] 0 A[So] 0 A[Sa,*.3])(0 o TTi o ((V'j,^)“\ Ca)), 
for all (f G Lf{QtA- 
The set B* is dual to B, that is 

^a{Ba>) = Saa', VA,A'gA. 
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Proof. To prove this theorem, we need to show that (33) is fulfilled. From the 


dual compatibility condition, it follows that (33) is fulfilled for all A G Ai and 


A' G A, where Wi is a structured chart. For unstructured charts Wi, it is obvious 
that 

= ^a^a'i VAi G Ai, and VA^ G Ai U (A\..4i). 

What remains to be shown is that = 0 for all A' G Ai\Ai. This 

becomes clear, as for all unstructured functions Ba the functional Aa evaluates 
to zero for all piecewise polynomials that are zero at the extraordinary features. 
This concludes the proof. □ 

In the following section we show approximation properties for a simplified 
configuration. Note that the results may be extended to more general configu¬ 
rations. 


5 . 4 . Approximation properties of spline manifolds on uniform meshes 

Having a stable dual basis, one can define a projection operator and use it 
to prove approximation properties of the B-spline manifold space. 

Here we consider h-refinement: a spline manifold H is given, with an initial 
mesh Tiio over it. From this we construct a family of meshes {Th} by mesh 
refining each structured chart. The refinement has to be consistent with the 
manifold structure, fulfilling Definition]^ in particular ([^, and Definitions [7j|^ 
Note that during h-refinement the structured charts and transition functions 
between them are unchanged, while the unstructured charts are modified in 
order to fulfill the mentioned Definitions. Accordingly, we assume that a family 
of nested spaces {Sh} is given. The subscript h always denotes the dependence 
on the refinement level. 

A physical manifold E is given on the initial mesh through a parametrization 


G G {Shor, as in (18). For simplicity, we consider a bivariate D, then E is a 


closed surface in the space During h-refinement G is kept unchanged. This 


gives a sequence of nested isogeometric spaces Vh via (19). 


To study the approximation properties of Vh we follow the approach in, e.g., 
El 11 HI] for structured spline spaces, which applies to the present framework 
as well. To keep it simple, we focus on a simple configuration: we consider 
a manifold obtained merging tensor-product patches with continuity but 

around the extraordinary vertices, where the continuity is only C^. This con¬ 
figuration is referred to as multi-patch B-splines with enhanced smoothness as 
in [lOj . where it has been first introduced and studied in the context of isogeo¬ 
metric analysis. On the coarsest mesh the length of the G° lines is p -I- 1 
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element edges, which corresponds to the condition that a function in {5/ip} is 
across the patch interfaces if and only if it has an extraordinary vertex in the 
closure of its support. The structured charts for this configuration can be taken 
as the union of each pair of patches that have a common interface. This is not 
the only possibility but it simplifies the next steps. In this case each transition 
function between structured charts is the composition of a translation and 
a rotation by a multiple of ^. To further simplify, we assume the mesh on each 
structured chart is uniform with mesh-size h. 

We assume by construction that each space Sh is complete. By this, we mean 
that for every structured chart uji the set Bi spans all piecewise polynomials of 
degree p and continuity stated above. 

The space is defined as 

L^(r2) = L^iuJi) for all structured charts Wi} , 

having the corresponding norm 

\\Hh(u)= Y. ( 37 ) 

QeTh 

where uJi Q is a structured chart covering Q, i.e. Q C tt^ Q(a;i g). Due to the 
isometry of the transition function, the L^-norm of a function defined on a chart 
fulfills || 7 ’|li 2 (^. ° i) structured charts Wi, Wj. Hence, 

the L^-norm on H is well-defined. Moreover, the definition of the L^-norm is 
independent of the level of refinement. We define bent Sobolev spaces 77^(D) 
in the same fashion. Bent Sobolev spaces are piecewise Sobolev spaces with 
some regularity at the element interfaces, see HHIIIZ]- For example, the space 
77^^^ (D) is defined as the closure of the space of piecewise C°° functions having 
the same continuity at the mesh lines of the space S, with respect to the norm 

p+i 

ll'^llwp+i(n) = 117^111,2(0) + Yh (3^) 

k=l 

where 

( 39 ) 

QeTh 

and I • is the usual fc-th order Sobolev seminorm. Obviously, all the spaces 

and norms can be defined accordingly on subdomains of D and are independent 
of the level of refinement. 

Using the dual basis defined in Theoremj^we introduce a projection operator 
: LF‘{VI) —>■ Sh onto the B-spline manifold space via 

</"-t = Y Aa((/')Ha (40) 

A.eAh 
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The projector is stable, uniformly with respect to h, i.e. 


(41) 

which follows directly from the stability of the dual basis as defined in The¬ 
orem see for example m- To prove its approximation properties we use the 
following Lemma. 

Lemma 3. Let Sr be the space of all piecewise polynomials with respect to a 
uniform Cartesian mesh of a box R GMf, such that 

• the mesh is formed by up to 2p + 1 elements per direction with meshsize 
h; 

• the polynomial degree is p in each direction; and 

• the continuity is globally with the exception of a mesh line e where 

the continuity is only i.e., Sr C C^~^{R \ e) n 

Let he the bent Sobolev space associated to Sr. Then for all (f> G 

there exists a a G Sr such that 

\\4> - (^\\l2{r) < ChP^^\(j)\Rp+i;R) (42) 

with a constant C only dependent on p. 

Proof. The size of R depends on the number of elements n = (ni,n 2 ) in each 
direction and the element size h. Given cf G TL^^^iR) there exists indeed a G Sr 
such that 

\\(t> - < C{p, h, n, e)\(j)\Rp+i;R) 

with C{p, h, n, e) independent of (p. The proof is the same as for the classical 
Bramble-Hilbert lemma, see, e.g., [55]. The dependence of the constant with 
respect to h, that is C{p,h,n,e) = C{p,n,e)hP~^^, follows from a scaling argu¬ 
ment. Finally, there are a finite number of different configurations for n and e, 
therefore we can set C{p) = max„,e C{p, n, e). □ 

See Figurej^for a possible configuration of R. Here, the line of (7° continuity 
is shown in blue. 

Given Q G Th, define Q C H in the following way: for each structured 
chart LOi, Q Ci Lli is the minimal box containing all the supports of functions 
in Bh whose support includes Q. We can now state the local approximation 
estimate. 
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Figure 5: Possible configuration of a box R and mesh line e for degree p = 2 

Theorem 4. Under the assumptions of this Section, for (p G it holds 

U-TIsM)\\lhq) < (43) 

where the constant C depends only on p. 

Proof. The first step of the proof is to show that, given (j) G there is 

an s S such that 

~ s \\ l 2 (^ q ) < ChP'^ l'/'l«p+i(Q)' (44) 

We are in one of two cases, either 

(a) Q contains an extraordinary vertex, or 

(b) Q C 7Ti{uJi) where uji is a structured chart. 




Figure 6: Different types of support extensions Q for p = 2 


See Figure for a representation of the two possible cases. In the figure the 
C'^-continuity lines are depicted in blue. The dark green element represents 
Q and the light green region represents its support extension Q. Due to the 
assumption on the length of the (7° lines, case 1) occurs when Q is adjacent 
to an extraordinary vertex, i.e., there exists an unstructured chart uJi such that 
Q C TTiioji). In this case we can split Q into Qi, ■ ■ ■ ,Qki such that each 


28 






























intersects only the segment 'Ki{si^i) (see DefinitionEach Su \q^ is a standard 
tensor-product spline space therefore we can use Ezi Section 2.2.2 and Section 
4.4] in order to construct splines G Sh \q^ which approximate (j) \q^ and match 
continuously in the whole Q. Since on the coarsest mesh the C® lines at the 
extraordinary vertex cover p-I- 1 element edges (see Figure]^, each interface 
between two adjacent sets of Qi,..., is covered by the continuity lines 
of s S and (441 follows. In Case 2), we can use Lemma then (441 follows 
by m. 


Having (44) and recalling the stability of the projector H^^, (43) is derived 
in the usual way, i.e. 


U-^sM)\\l^q) = U - (T - ns^{(j) - a)\\L2(Q) 


The stability constant as well as the approximation constant only depend on p, 
which concludes the proof. □ 


By means of inverse estimates and generalising (441, (43) can be extended 
to higher order Sobolev norms. For 0 < q < p + 1 it holds 


\cl,-IlsM)\H.iQ) < (45) 

where the constant C depends only on p and q. The details are not reported 
for the sake of brevity. 

We can extend the result to a mapped domain E, in this case a closed surface 
in parametrized by G e {Shg)^- We assume that G is regular, that is, there 
exist constants c, c, with 


c > det(VG^(x)VG(x)) > c > 0 


for all X S Q and for all Q G Tho- Note that, in general, we cannot define 
Sobolev spaces of any order on E, due to the lack of smoothness of the manifold 
itself. However the space on E can be defined as 

L^{^) = {f:foGGL^n)} 


and the corresponding norm is given via 

II/IIl^(e) = ||/oG(det(VG^VG))i/4|U2(a). 


29 









The bent Sobolev spaces can be defined similarly, as in (38)-(39). See [301 in] 
for more details about Sobolev spaces on manifolds. Then, for all / G Lp‘{Y?}, 
we can define the isogeometric projector 


nv.(/) = n5,(/oG)oG 


-1 


Approximation properties of easily follow from the ones of stated in 
Theorem following the same approach of [551 131 ] • 

Theorem 5. Under the assumptions of this Section, for all f G 

11/ - nvJ/)IU2(s) < (46) 

where the constant C depends only on G and p. 

Note that all the results presented here extend naturally to volumetric do¬ 
mains. In that case, the lines of G° continuity extend to faces of G° continuity 
in a vicinity of the extraordinary vertices and edges. 


6. Conclusion and possible extensions 

We have introduced a new general mathematical framework, based on man¬ 
ifolds, for the definition and the analysis of unstructured spline spaces. As it 
is done in ^ the main idea is to decompose the domain into charts, which 
are meshed with quadrilaterals or hexahedra. Then spline basis functions and 
dual functionals can be defined locally on each chart. Unstructured charts are 
necessary to cover extraordinary vertices and edges of the domain. 

We have used this framework to generalize the dual-compatibility condition 
of [23112] to unstructured spline spaces, and in particular to analyze the ap¬ 
proximation properties of splines with high continuity everywhere except in the 
vicinity of extraordinary vertices and edges, where the continuity is only . 
Although the analysis was restricted to the low continuity case, the framework 
allows for the definition of spline functions with higher smoothness, and their 
analysis will be the aim of future work. 

In our definitions the physical domain is necessarily a manifold. However, 
since we are defining the charts in the parametric domain, and not in the physical 
domain, it is possible to extend our framework to non-manifold domains using 
special bifurcation charts (such as T-shaped or X-shaped charts for curves, etc.) 
This could be of interest for certain beam or shell formulations, or for the proper 
representation of the medial axis or medial surface of an object, for instance. 

Finally, for the sake of simplicity we have restricted ourselves to B-splines 
and T-splines on quadrilateral/hexahedral meshes. The framework can be easily 
generalized to other spline spaces, such as NURBS or trigonometric splines. 
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Appendix A. Spline manifolds with boundary 

We can extend the definition of a spline manifold to a manifold with bound¬ 
ary. To do so, we first need to extend Definition defining a suitable proto¬ 
manifold that takes into account the boundary. 

Definition 24 (Manifolds with boundary). A proto-manifold with boundary is 
a generalization of a proto-manifold, where the charts {wi}i=i,...,Ar are given 
as uji = interior(iOi) U 7 ,^, such that interior(uji) are open polytopes forming 
a standard proto-manifold (with transition domains interior{ojij) and transi¬ 
tion functions ipij) and each 7 ^ C duii is a part of the boundary of the chart 
interior{(jJi). Moreover, the transition domains fulfill ujij = interior 
with 7 i j = "findujij, and the transition functions are the continuous extensions 
of tpij and map 7^7 onto 7^7 and interior (uiij) onto interior 

Similar to the standard parameter manifold in Definition we can define 
the parameter manifold with boundary Q via the equivalence relation induced by 
the transition functions. Since the transition functions always map the interior 
onto the interior and the boundary onto the boundary, the parameter manifold 
D can be separated into interior {it) and the boundary denoted by P. 

This definition of the boundary of the (open) parameter manifold is equiv¬ 
alent to the classical definition of a manifold with boundary, as discussed in 
[T]. In this case every boundary point of the manifold has a neighborhood that 
is homeomorphic to the half d-ball. Note that the boundary P itself can be 
interpreted as a topological manifold of dimension d — 1. 

For the definition of the spline spaces, we assume the following. 

Assumption 25. The local boundary 7 ^ is conforming with respect to the 
elements, i.e. there exists a subset of faces of elements q G Ti that forms a mesh 
for the boundary 7 ^. 

To be able to define manifolds with boundary containing non-convex fea¬ 
tures, we need additional types of charts, so called boundary charts. To avoid 
the tedious formal definition of boundary charts, we present figures that should 
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explain the ideas behind. In Figure 1^ we present different types of two- 
dimensional boundary vertices, regular boundary vertices (a) - (c), as well as 
a non-regular boudary vertex (d). The vertices in Figures [7(a)| and |7(b)| can 
be covered by the boundaries of structured charts (see Figure [8(a)] - |8(c)] ). For 
the vertex in Figure 7(c) we need a special boundary chart (see Figure [8(^ , 
associated with the function which is non-zero at the corner. Note that the 
boundary vertex in Figure [7(d)| is discarded since it does not allow for functions 
that are non-zero at the corner. 


(a) (b) (c) (d) 


Figure A.7: Different types of valid (a)-(c) and non-valid (d) boundary vertices in 2D 


(a) (b) (c) (d) 


Figure A.8: Structured charts with boundary (a)-(c) and boundary chart (d) in 2D 


In the three-dimensional case, there are more different configurations to con¬ 
sider, which are listed in Table A.3 Here, we need to consider two different types 
of boundary charts. 


structured edge 
structured vertex 

unstructured vertex 


regular boundary edge 

regular boundary vertex 

hanging boundary vertex 

partially unstructured boundary vertex 


Table A.3: Classification of valid boundary edges and vertices in 3D 


Figures [A.9| and |A.10| depict several possible three-dimensional charts with 
boundary. In Figures |9(a)|j9(c) we show several examples of structured charts 
with boundary. Figure [9(d)| shows an unstructured edge chart with boundary. 
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In Figure |A.10 we show the two different types of boundary charts, that are 


9 8901 


(a) 


(b) 


(c) 


(d) 


Figure A.9: Structured charts with boundary (a)-(c) and unstructured edge chart with bound¬ 
ary (d) in 3D 

needed to represent all meshes of practical interest. The one in Figure |10(a)| 
is a tensor-product of a two-dimensional boundary chart with an interval in 
the third direction. The chart depicted in Figure |10(b)| is a boundary chart 
corresponding to a non-convex vertex at the boundary. One can include more 
complex boundary configurations by introducing unstructured boundary charts, 
which we will not consider for the sake of simplicity. 




(a) (b) 


Figure A. 10: Boundary charts in 3D 


Concerning spline manifold spaces on parameter manifolds with boundary, 
we need to adjust the condition on proto-basis functions in equation Q to 
obtain 

lim 6 a,(C)= 0. (A.l) 

With this modification, the spline manifold space can interpolate at the bound¬ 
ary. The theory concerning dual-compatibility and approximation properties 
presented in Section can be generalized directly to spline manifolds with 
boundary. Most importantly, both Theorems and extend directly to mani¬ 
folds with boundary, using the boundary charts introduced here. Thus extending 
the approximation error bounds to manifolds with boundary of general topology. 
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